Dynamical models for sand ripples beneath surface waves 
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We introduce order parameter models for describing the dynamics of sand ripple patterns under 
oscillatory flow. A crucial ingredient of these models is the mass transport between adjacent ripples, 
which we obtain from detailed numerical simulations for a range of ripple sizes. Using this mass 
transport function, our models predict the existence of a stable band of wavenumbers limited by 
secondary instabilities. Small ripples coarsen in our models and this process leads to a sharply 
selected final wavenumber, in agreement with experimental observations. 



I. INTRODUCTION 

When a flat surface of sand is exposed to the flow of 
air or water, patterns known as ripples, dunes, sandwaves 
and draas are formed |fl-p|. Here we focus on the so- 
called vortex ripples jp (Fig. 0) which are created by os- 
cillatory fluid flow, e.g., beneath surface waves. Ripples 
are of interest to coastal engineers since their properties 
determine the friction of the flow in the coastal region 
jj] |To| , the dissipation of surface waves |ll]] and the net 
sediment transport over the ripples |12|. More recently 
ripples have attracted the attention of physicists inter- 
ested in non-equilibrium systems |^|||L3|-^0) . 

The physics underlying sand ripple formation involves 
the interaction between the turbulent fluid flow and a 
granular medium, and is therefore extremely complex. 
A description of the pattern forming aspects is hindered 
by the strong nonlinearity of the fully developed ripples 
due to the subcritical nature of the initial bifurcation 
from a flat bed. Previous theoretical studies of this ini- 
tial bifurcation ]2l|-|2^| have described the onset of ripple 
formation. Vortex ripple pattern formation occurs, how- 
ever, far from equilibrium: typical wavelengths of fully 
developed ripples can be a factor five larger than those 
predicted by (weakly nonlinear) analysis p3). 

In this paper we will discuss the pattern forming as- 
pects of fully developed vortex ripples. Many of the prob- 
lems associated with the complicated underlying phe- 
nomenology can be circumvented by noting that the sizes 
of the ripples are the most relevant parameter for de- 
termination of their dynamics; further details of their 
shapes are not important. Dynamical equations for the 
evolution of the ripples can then be constructed once the 
mass exchange between ripples of certain sizes is known. 
We base our expression for this mass exchange on de- 
tailed numerical simulations of the flow and sand trans- 
port over vortex ripples (see below), hence going beyond 
a pure "toy-model" approach. As far as we are aware, 
the model presented here is the first to capture both 



FIG. 1. Side-view of a vortex ripple pattern under oscilla- 
tory flow in a long, slender channel (courtesy of J.L. Hansen). 
The arrow indicates the amplitude of the fluid oscillations in 
the bulk. 

instabilities and coarsening of fully developed vortex rip- 
ples. 

The outline of the paper is as follows. We start with 
a brief description of the main phenomenology of ripples 
in section 0. Although the amplitude of the fluid os- 
cillations determines th e len gth of the ripples, a dimen- 
sional analysis (section II A ) reveals that the most rele- 
vant dimensionless control parameter is the Shields pa- 
rameter which characterizes stress at the sandy surface. 
We discuss our numerical simulations of the mass ex- 



change between vortex ripples in section II B . Section III 
is devoted to the formulation of simple ripple models in 
one-dimensional geometries. The linear stability of these 
models is performed in section [II B , and the coarsening 



and selection of the final ripple patterns starting from 
random initial conditions is discussed in section HID. In 
section |y| we extend our model to two dimensions and 
discuss the impact of defect motion on the selection of 
the final two-dimensional pattern. 



II. VORTEX RIPPLES 



Following the much earlier work of Ayrton p5[ , the 
study of vortex ripples was taken up again by Bagnold 
in 1946 0]. In this seminal study, Bagnold distinguished 
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between rolling grain ripples and vortex ripples. The for- 
mer are generated when starting from an unstable flat 
bed [[| and consist of small triangular ridges separated 
by a comparatively long stretch of flat bed. These rolling 
grain ripples grow and coarsen to become vortex ripples 
with no flat bed between them. Here the flow is dom- 
inated by separation bubbles (vortices) on the lee sides 
of the fully developed ripples. We will concentrate on 
these fully developed vortex ripples, since recent studies 
have confirmed Jl^,^(| that rolling grain ripples essen- 
tially constitute a transient. 

Many experiments have studied the average wave 
length of fully developed ripples p6|-|2"i| as a function 
of, e.g., the amplitude and frequency of the fluid motion. 
It appears that the (dimensional) length of the ripples 
(Adim) is proportional to the amplitude of the oscillatory 
flow (a), and roughly independent of its frequency. Es- 
timates in the literature of the proportionality constant 
^dim/o- range from one to two, with a preference for val- 
ues around 1.3 (see Fig. 8 in p7[|). 

Recently the ripples have also been studied from the 
view point of pattern formation. Both Scherer et al. |2^] 
and Stegner and Wesfreid |ll| studied a one-dimensional 
annular system in which the conservation of sand is guar- 
anteed. Stegner and Wesfried |l9| observed strong hys- 
teresis when the driving amplitude of fully developed rip- 
ples was ramped up and down: an increase in the am- 
plitude of the driving yielded larger ripples, while for a 
decrease, the ripples did not change length. Lofquist also 
observed hysteretic behavior, but in this case the rip- 
ples were initially stable for both an increase or decrease 
of the driving amplitude (t]] . Hysteresis of the ripples was 
also observed in a recent set of field measurements l29| . 



A. Dimensional analysis and setup of the problem 

Ripples are governed by a large number of dimensional 
parameters which characterize fluid flow and sand. We 
will show that while in general three dimensionless pa- 
rameters (density ratio of fluid and sand grains s, set- 
tling velocity w s and maximum Shields parameter 9 max) 
characterize the system, for the case of interest here 
(sand/water systems in the regime where suspension is 
unimportant) the only free parameter is the Shields pa- 
rameter. 

Ripple formation is driven by an oscillatory fluid mo- 
tion with amplitude a and angular frequency lo. The 
Reynolds number Re for this situation is a 2 aj/v, where 
v is the fluid viscosity. For water in a typical experi- 
mental situation (a = 5 cm, u> — 3 s" 1 ), Re is of order 
10 3 , and the flow is turbulent. Therefore, large scale flow 
structures such as separation bubbles are independent of 
the Reynolds number and hence viscosity. For turbulent 
flow, the roughness of the bed is of minor importance as 
long as the typical grain sizes are much smaller than a. 
The only relevant length scale is then a, which we use to 
define the non-dimensional ripple length as A — Xdim/a. 
The large scale flow is then completely specified by the 
boundary conditions, i.e., the shape of the ripples. 

The sand introduces four new dimensional parameters 
into the problem. These are, respectively, the density of 
water p w and sand p s , the median diameter of the grains 
d and gravity g. From these we form the following three 
non-dimensional parameters: 
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FIG. 2. The flow over a ripple at three instants in the first half of the wave period. The system consists of a single ripple 
in a system with periodic boundary conditions and length A = 1.15. The bottom three figures show the spatial profile of the 
Shields parameter on the bed at the corresponding times. For this C&SG Umax — 0.13. 
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The relative density of the grains s has a value of 2.65 
for quartz sand in water. The settling velocity lu s char- 
acterizes the amount of sand kept in suspension; here 
we assume a regime where the settling velocity is large 
(w s >0.15) such that suspension is not important. This 
leaves us with the last parameter, 6, 
which is known as the Shields parameter. The Shields pa- 
rameter expresses the ratio between the drag and gravita- 
tional forces on a single grain and depends on the shear 
stress Tbed{x,t), which varies with time and along the 
profile of the ripple. Following |3(| we propose to use 
the maximum shear stress on a flat bed, T max , to char- 
acterize the flow. For laminar flow T max can be found 
exactly from the solution of Stokes' second problem |sT| . 
For turbulent conditions, which prevail here, an analyti- 
cal expression does not exist. We will follow coastal en- 
gineers in using an empirical relation for the maximum 
shear stress 1321 : 
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(2) 



Note that the instantaneous Shields parameter on a rip- 
pled bed 9(x, t) can be several times larger than 9 max . 

The transport of sand q takes place in a thin layer 
above the bed, the so-called bed load layer (for an intro- 
duction to sediment transport see chapter 7 i n ref. 
The non-dimensionalized flux of sand <\> = qj -J g(s — l)d 3 
in the bed load layer is a function of the local Shields 
parameter and can be modeled as: 



cj> = a{6- 9 C ) 



(3) 



When 8(x, t) smaller than a critical value 8 C for all x, 
which for turbulent boundary layers is approximately 
0.06 |32[ , sand grains do not move and the ripple pro- 
file freezes. The constants a and [3 have been determined 
empirically by Meyer-Peter and Miiller Q to be approx- 
imately a = 8 and /5 = 1.5, which are in good agreement 
with theoretical estimates |54|,[55]]. The formation and 
the dynamical properties of the ripples are mainly deter- 
mined by the fluid flow, so the exact values of the con- 
stants a and /3 together with the detailed form of Eq. (||), 
turn out to be relatively unimportant for the content of 
this work. 



B. Numerical studies and mass transport 

The computational model that we have developed to 
study the ripples calculates turbulent fluid flow over rip- 
ples based on the standard k-ui turbulence model |5(|[37| ■ 
Once this flow is known, the sediment transport, which is 
governed by the shear stress on the bed, can be calculated 
from Eq. mh. 

In Fig. | we show some results for the flow and the 
non-dimensionalized shear stress 4> for A = 1.15. We see 
that 




FIG. 3. The spatial profile of the net mass transport in 
the first half of the wave period (/J 1 ^" 4>(x,t)dt) for ripples of 
length A = 0.6, 1.0 and 2.0. Note that the ripple trough is 
located in the middle of the figure. (9 max = 0.13). 

there are two mechanisms that generate the shear stress 
on the bed, namely the converging flow on the "wind" 
side of the ripple (here left) and the separation bubble 
formed in the lee (right) side. Typically, these stresses 
are several times stronger than the stresses on a flat bed. 
The separation bubble, where the flow near the bed is 
directed opposite to the mean flow direction, is clearly 
visible. This bubble moves out into the trough of the 
ripple (ujt — 90°), where it stays (ujt — 150°) until it is 
thrown over the crest as the flow reverses. 

The shear stresses are uphill on both sides of the ripple, 
and consequently sand is transported from the trough 
toward the crest; the result is a steepening of the ripple 
profile. This steepening continues until the slopes of the 
ripple reach the angle of repose, when avalanches limit 
the growth of the ripple slopes. As a consequence, most 
slopes of the fully developed ripples are close to the angle 
of repose. These fully developed ripples are thus approx- 
imately triangular, joined by smooth troughs, which is 
also evident from experiments p9| . 

Ripples interact by exchanging sand between their 
neighbours over the troughs. The amount of this mass 
flow is closely connected to the extension and strength of 
the separation bubble. We have studied this mass trans- 
port as a function of the non-dimensional ripple size A. In 
Fig. [3] we show the net sediment transport during the first 
half wave period, for short (A = 0.6), medium (A = 1.0) 
and large ripples (A = 2.0). For short ripples the sepa- 
ration bubble almost covers the space between the two 
ripple crests, but it is not very strong, giving rise to a 
small transport. For long ripples, the separation bubble 
does not reach over the trough, again giving only a small 
mass exchange between adjacent ripples. Most mass is 
exchanged for medium sized ripples, where the separation 
bubble is both strong and reaches over the trough. 

We define / as the amount of sand transported over 
the trough during the first half wave period: 
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FIG. 4. /(A) for 9 max = 0.075, 0.15, 0.23, 0.30, 0.45, 0.60 
and 0.75, scaled with /(A = 1.0). The inset show 10]^ ax (full 
line) and /(A = 1.0) (dashed line) in double log scaling. 



/(A) 



7T / ~U 



(x tr ,t)dt, 



(4) 



where Xt r is the position of the trough. The minus sign 
is simply related to the fact that the fluid and mass flows 
have opposite directions during each half period; here we 
wish to have a positive /(A). 

The rescaled mass exchange /(A)//(A = 1.0) is shown 
in Fig. U for values of 9 max ranging from 0.075 to 0.75. 
The rescaled graphs of the mass exchange (i) collapse 
in good approximation, and (it) have a single maximum 
around A = 1.0. In our model, developed in section III 
below, we will incorporate these two properties. If the 
critical Shields parameter had been zero, the rescaling 
factor /(A = 1.0) would have been proportional to fl^as- 
That this is almost, but not exactly, the case is seen in 
the inset in Fig. 0. 



III. DISCRETE MODELS FOR 
ONE-DIMENSIONAL RIPPLES 



In this section we will introduce and study simple mod- 
els for ripples in one-dimensional geometries. We assume 
that the angles of the ripple slopes are fixed, so that the 
only degrees of freedom are the lengths of their left and 
right slopes. Two different versions of the model will be 
described. In the simplest case we only take the total 
ripple sizes Ai into account (see Fig. ||). The ensuing 
"minimal model" is formulated in section III A. and is 



analyzed theoretically in section HI B . A more refined 
model which takes the lengths of le ft and right slopes 
into account is presented in section III C (see Fig. g), 
and nu merica l simulations of this model are presented in 
section HID. 
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FIG. 5. An example of how the transport in the trough be- 
tween two ripples is constructed as the transport in the two 
half periods. 



A. Minimal model 

In this model ripples are triangular and symmetric and 
characterized by their length A^ . We will now determine 
the mass transfer between two ripples with lengths Ai 
and A2 (Fig. ||) from the information that we have for the 
mass transfer between equal ripples. When Ai and A2 are 
approximately equal, one expects the size and strength 
of the separation bubble emanating from the crest of rip- 
ple 1 to be independent of the size of ripple 2. This is 
our central assumption: the mass transport during a half 
period only depends on the size of the ripple that creates 
the separation bubble. 

Let us denote the first half period of the driving, when 
the flow is from left to right, by a subscript J, and the 
second half by II. Under our assumption stated above, 
we obtain: Amj = f(Xi) and Amu = — /(A2), where 
Amj denotes the change in the mass of ripple 1 in the 
first half period. During each half period, the amount of 
mass transported is small in comparison to the mass of 
a single ripple. We therefore can neglect changes in the 
ripple shapes during a half period, and obtain the mass 
flow during a full period, Am, by simply adding up the 
half period mass-flows: 



Am = /(Ai)-/(A a ) 



(5) 



Clearly Eq. (§) can be extended to the case of a row 
of ripples. Then the mass- flow to ripple i, Arm, is 
due to interactions with both ripple i — 1 and i + 1: 
Am 1 = 2/(A,)-/(A t+1 )-/(V 1 ). 

To close the equations we need to relate the mass flow 
to a change in the size of the ripples. Since the mass-flow 
is small, it is reasonable to assume that the change in 
ripple size is linear in the mass transport. The greatest 
simplification is obtained if we assume all ripples to be of 
near equal size, so that the ratio Am/ AX is equal for all 
ripples. Taking the continuum time limit and rescaling 
time to absorb a proportionality constant we obtain: 



dXi/dt = -f(Xi-i) + 2f{Xi) - /(A i+X ) 



(6) 



4 



The total length of a system of ripples evolving according 
to Eq. M is conserved, but the total mass is not; we will 



U 



discuss this further in section III C 



Finally, we supplement the model with an annihilation 
rule which removes ripples that have shrunk to size zero, 
and a creation rule, which adds ripples in the troughs be- 
tween ripples of sizes larger than a certain length X max 
that will be specified in the next section. 



B. Equilibria and stability 

There are three types of equilibria in the minimal 
model (o). (i) Homogeneous states where all A's are 
equal, (u) "Period two" states, for which /(Aj) = /(Aj+i) 
but \ ^ Aj+i (see Fig. 16 in |2(| for similar states). 
(Hi) More complicated equilibria constructed by arbi- 
trary juxtapositions of ripples of lengths A a or A& when 
/(A a )=/(A fe ). 

The linear stability of the homogeneous state follows 



from setting A^ = X eq + Si and linearizing Eq. 



dS i /dt = -f'(\ eq )(6 i -i 



25,. 



(7) 



This is the linear stability equation for the space- 
discrctizcd diffusion equation, with diffusion coefficient 
— f'(X eq ), and the sign of f will be important. As we 
demonstrated in section II B , the mass transport / dis- 
plays a single maximum as a function of A at a value that 
we will refer to as \ m in- When X eq is larger (smaller) 
than X m in, —f'(Xeq) is positive (negative) and the pat- 
tern is stable (unstable). Hence the smallest possible sta- 
ble wavelength is \ m im where / has a maximum. This 
instability can be seen directly from the mass transport: 
when we inspect two unequal adjacent ripples with sizes 
larger than X m i n we obtain from Eq. (^|) that mass will 
flow from the larger to the smaller ripple, hence leading 
to a stable equilibrium, while if their sizes are smaller 
than X rn in mass flows from the smaller to the larger rip- 
ple, leading to an instability. 

An additional instability occurs for large ripples when 
their troughs lie outside the separation zone (see A = 2.0 
in Fig. in this case the flow creates new ripples in 
the troughs. This instability has been observed in ex- 
periments and also in our numerical studies f37| . 
This instability is consistent with our models when we 
assume that / is defined for arbitrary small ripples. For 
a homogeneous pattern of large ripples where f{X eq ) < 
/(0), infinitesimal ripples inserted between the large rip- 
ples will gain mass and grow, and the maximum value 
Xmax where homogeneous patterns are stable, is given 
by f(X rnax ) — f(0). This is the motivation for having a 
creation rule in the model. 

The period-two and more complicated equilibrium 
states can be shown to be unstable in our framework |38| . 
Thus our model illustrates an important consequence of 
the shape of the mass exchange function. There is a 
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FIG. 6. A sketch of the ripple profile with the triangular 
ripples and the quantities used to describe the ripples. Note 
that the ripples do not need to be symmetric. Also shown 
is the exchange of sediment from the right ripple to the left 
ripple during the first half wave period. 

band of wavelengths for which ripple patterns are sta- 
ble; outside this band, short wavelength instabilities oc- 



C. Refined model 

Both our numerical studies and experiments [fl9f fre- 
quently display ripples that are asymmetric during their 
evolution (although, on average, they are not). We ex- 
tend the minimal model from the previous section to al- 
low for asymmetric ripples by characterizing the ripples 
by the length of both their left (L) and right (r^) slopes; 
obviously A.; = U + (see Fig. ||). In addition, such a 
model can be tuned so as to conserve mass. 

As before we assume that the lee side of ripples de- 
termines the size of the separation bubbles. During the 
first half period the bubble takes Am/ mass from the left 
slope of ripple i + 1, and transports this mass to ripple 
i; the ratio between the mass deposited on the left and 
right slopes of ripple i is given by a parameter a that we 
always fix at a value of 0.5 (see Fig. ||). The mass flow 
in the first half period is therefore: 



Am,,,/ = -/(2ri_i) + (1 - a)f(2n) 
Am r ij = o-f(2n) . 



(8) 



The mass-flow in the second half period follows by sym- 
metry. Assuming that the mass transport is small, we 
can neglect the change in ripple size during one half cy- 
cle, and add the contributions from each half period. 

To obtain a closed set of dynamical equations we have 
to establish how the lengths k and evolve under a cer- 
tain mass flow. When an amount of mass M is deposited 
on the right slope of ripple i, we incorporate this by an in- 
crease of li and a subsequent decrease of U+i] the length 

itself does not change. Assuming for simplicity the 
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FIG. 7. The interaction functions used in the models. 
Smooth function, (full line) and the bi-linear function (broken 
line), for a value of A ma:E = 1.5. 

angle of repose to be 45°, we can calculate the volume 
of the slab of deposited sand and find that the change in 
the length of nearby ripples is: 

Ak = Am„/(2r t ) (9) 
Al l+1 = -Am ri /(2r i ) . (10) 

Ignoring higher order effects we obtain the total change 
in the length as a function of the mass changes as: 

Am H+1 , Amu 



An = - 
AL = 



2k 
Am r i 



2n 



-l 2Z j 
Am r ,;_i 

~ 2n-i 



(11) 



This relation together with the mass flow from (|J) de- 
fines the refined model. This model has the same linear 
stability properties as the minimal model defined in Eq. 

The total length of the system is conserved, and the to- 
tal mass is approximately conserved. The masses which 
are ignored are associated with the small areas that are 
cross-hatched in Fig. |(| It is possible to formulate the 
model in a strictly mass-conserving manner, by updat- 
ing the slope lengths when both removing and depositing 
mass, but this does not alter the model in any substantial 
way. 

In our numerical simulations two different forms of the 
mass transport function /(A) were used. Both functions 
have a maximum at A = A. m ; n = 1/2 and are zero at 
and X m ax- The simplest function that satisfies these 
requirements is bi-linear, while a smooth function with 
a quadratic maximum that satisfies /(0) = /(A max ) = 
can be constructed as the sum of a linear function and a 
square root (see Fig. fj]): 
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FIG. 8. An example of the dynamics of the model using the 
linear interaction function and A m aa; = l-35 The time scale is 
arbitrary. 



1.6 resembles the one 



This smooth function for A m , 
found from the computational flow model in section II B 



D. Coarsening of fully developed ripples 



When ripples are grown experimentally from a flat bed, 
initially many small ripples are created. They subse- 
quently coarsen and form a final regular steady state with 
a well-defined final wavelength (see for example Fig. 1 
in Our model shows the same behavior for initial 

conditions of (disordered) unstable small wavelength pat- 
terns. An example of such evolution is shown in Fig. [§. 
A fast coarsening process is seen in the beginning (t < 1), 
followed by a slower relaxation toward an equilibrium 
state. The important dynamical process leading to the 
equilibrium state is the annihilation of ripples, with each 
annihilation resulting in a longer average ripple length; 
creation does not play a role here. After the final anni- 
hilation, slow diffusive dynamics sets in. 

The stability analysis performed in section 
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shows 

that a wide range of ripple wave lengths can be linearly 
stable, namely A TO j„ < A < X m ax- We will show here 
that, starting from small ripples, the dynamics leads to 
the selection of a sharply defined final wavelength. We 
assume periodic boundary conditions in our simulations. 
The parameters entering the model are the length of the 
domain L and the maximum ripple length X m ax- The 
initial conditions are disordered ripples with an average 
wavelength Ao<A m j n . 

We found that the final wavelength is quite indepen- 
dent upon the initial average wavelength Ao (when this is 
sufficiently small) and the initial degree of disorder. This 
result could not have been predicted a priori from the 
model equations, but is in good agreement with experi- 
mental evidence. 
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FIG. 9. The equilibrium wavelength as a function of X, nax . 
The squares correspond to results obtained using the bi-linear 
interaction function, the open circles corresponds to the 
smooth function and the filled circles corresponds to the 2D 
model. The dashed line indicates the maximum possible 
wavelength, A = Xmax- The initial number of ripples were 
1200, and the initial ripple length was 0.7±0.05. 

The final wavelength does, however, depend on the shape 
of the interaction function and the value of X max - In 
Fig. ^, X eq is plotted as a function of X max for the two 
interaction functions. The final wavelength appears to be 
a nontrivial function of X ma x for both interaction func- 
tions. The interaction function which resembles the one 
from the numerical flow calculations (the smooth func- 
tion with X max = 1.6) results in an equilibrium wave- 
length of X eq = 1.28 ± 0.03, a result which is in good 
agreement with ripple lengths measured in experiments. 



IV. TWO-DIMENSIONAL RIPPLE PATTERNS 



Unless one forces the ripple patterns to occur in a 
narrow channel or annulus, ripple patterns are two- 
dimensional, even though the alignment perpendicular 
to the flow yields quasi one-dimensional patterns. How- 
ever, during the evolution toward the final state the pat- 
tern contains many defects 0]. We have thus extended 
our model to study their role in the selection of the final 
state. 

In our two-dimensional model the individual rows of 
ripples are similar to those in Fig. || and are labeled by 
an indices i and j, where j is the new coordinate perpen- 
dicular to the driving direction. In the z-direction, the 
mass flow is given by expression (|^) . We then determine 
which ripples are neighbors in the j-direction, and im- 
pose an angle of repose in the j direction by inducing a 
flow of mass between ripple j and j + 1 when their height 
difference is above a certain maximum. At a defect, such 
mass flow can nucleate a new ripple in the trough of an 
adjacent row. Finally, it is reasonable to expect that a 
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FIG. 10. An example of a climbing defect in a system with 
30 ripples in the middle and 31 ripples at the top and bottom. 
(Ai= 1.29, A 2 = 1.33). 

j-flow is induced when the ripples are not aligned perpen- 
dicular to the main flow, and the simplest choice for such 
flow between ripple j and j+1 is —Amy = C x {xj+i —Xj). 
The coupling in the j-direction is thus diffusive, and ba- 
sically acts to align the ripple crests perpendicular to 
the oscillatory motion. In the simulations which are pre- 
sented, the value of C x has been fixed to 0.08. The 
qualitative results are not sensitive to a change of this 
parameters. 

To study the motion of defects in this model, we ini- 
tiate the system with two patches of nearby wavelengths 
Ai and A 2 separated by two defects. The motion of 
these defects depends on the values of the wavelengths 
(Fig. [h]), and we find that when these are larger than 
1.28 ±0.02 the defect climbs in the direction of the low- 
est wavelength, otherwise it moves in the direction of the 
largest wavelength. Thus in a pattern with many de- 
fects, as encountered during the coarsening process, one 
expects that only regions with wavelengths larger than 
Xdef — 1-28 will survive. We can therefore expect the 
final wavelength X eq in a two-dimensional system to lie 
between X de f and X max - 

To check this we have performed simulations in large 
systems with initial conditions consisting of unstable 
ripple patterns of wavelength 0.5. The system rapidly 
coarsens and evolves to a state where most of the wave- 
lengths lie in the ID stable regime A m j„ < A < X max . 
The dynamics then slows down dramatically and be- 
comes dominated by defect climbing. In the final relaxed 
state of the system the peak of the distribution of ripples 
lengths lie between Xdef and X max . Runs with differ- 
ent values of X max have confirmed that the equilibrium 
wavelength in the two-dimensional model is systemati- 
cally larger than in the one-dimensional case (see Fig. |]). 
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V. DISCUSSION AND CONCLUSION 

By focusing on a realistic law for the mass exchange 
between adjacent ripples, simple models have been for- 
mulated that capture a number of phenomena observed 
in real ripples. First of all, our model predicts the ex- 
istence of a finite band of stable ripple wavenumbers 



A„ 



< A < A„ 



which is consistent with the hystere- 



sis observed in experiments [|7|]29|1 . The model predicts 
that the instabilities encountered once these boundaries 
are crossed are of short wavelength nature, in agreement 
with experiments on one- and two-dimensional sand pat- 
terns performed in Copenhagen 

Coarsening which occurs in the intermediate stages of 
the development of vortex ripples is present in our mod- 
els, and we predict that even though there is a finite band 
of stable ripple patterns, the dynamics selects a well de- 
fined final wave length. The exact value of the final wave 
length depends on the details of the mass exchange func- 
tion, however for a function similar to the results from 
our simulations of the fluid and sand flow, we find an 
equilibrium wavelength of A = 1.28. 

The good collapse of the mass exchange function with 
the maximum shear stress 9 max indicates that the final 
wavenumber should be independent of 6 max as long as 
suspension is not important. 

Following the picture of the mass exchange as in the 
models, it is clear that the maximum value of the mass 
exchange function sets a time scale for the evolution of 
the ripples. We have shown that this maximum value is 
approximately proportional to 6 max . Thus the time scale 
of the evolution of the ripples can be expected to scale 

as "max- 

Finally we have demonstrated that the models can be 
applied to two dimensional ripples patterns, and have 
found that defect motion renders the final wavelength of 
ripples in two-dimensions larger than in one dimension. 

All these predictions are open to experimental verifi- 
cation. In particular, we are eager to see how consistent 
the mass exchange mechanism proposed here is for real 
data of ripple evolution. 
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